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Abstract 

Magnetic properties of quark matter at finite temperature are discussed by evaluating the mag- 
netic susceptibility. Combining the microscopic calculation of the self-energy for quarks as well 
as the screening effects for gluons with Fermi-liquid theory in a consistent way, we study the 
temperature dependence of the magnetic susceptibility. The longitudinal gluons have the static 
screening given by the Debye mass, and have a standard temperature dependence of 0{T^). An 
anomalous In T term arises in the magnetic susceptibility as a novel non-Fermi-liquid effect 
due to the anomalous self-energy for quarks given by the dynamic screening for transverse glu- 
ons. We then extract the critical (Curie) temperature and present the magnetic phase diagram on 
the density-temperature plane. 
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1. Introduction 

Nowadays the phase diagram of QCD has been extensively studied theoretically by lattice 
calculations [1] and the effective models or experimentally by high-energy heavy-ion collisions 
(RHIC,LHC) |2] and observation of compact stars [3]. In particular, matter at high density but 
not high temperature should be interesting, since the Fermi surface is a good concept there and 
we can study interesting particle correlations such as superconducting pairing or instability of 
the Fermi surface in a clear way. It is also related with a quest of "new form of matter" inside 
compact stars or during their thermal evolution |4]. 

We discuss here magnetic properties of QCD by the use of Fermi-liquid theory IH S 7]. 
A possibility of spontaneous spin polarization has been suggested by one of the authors isHSfl, 
using the one-gluon-exchange (OGE) interaction. Higher-order calculations for the free energy 
have also supported this idea . If it is realized, it may give a microscopic origin of strong 

magnetic field observed in compact stars, especially magnetars 1 11]. The mechanism of sponta- 
neous magnetization is very much similar to the one for the electron gas 1 12, 13, 14], where the 
leading-order contribution comes from the Fock exchange interaction. Electrons with the same 
spin can effectively avoid the repulsive Coulomb interaction due to the Pauli principle. On the 
other hand, the kinetic energy is increased as the number of each spin is different. Hence spon- 
taneous polarization may occur at a peculiar density. Differences are color and flavor degrees of 
freedom for quark matter. 
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We apply the Landau Fermi-liquid theory (FLT) to elucidate the critical behavior of the mag- 
netic phase transition at finite density and temperature. We evaluate the magnetization of quark 
matter by applying a small magnetic field B. Then the magnetic susceptibility appears as a co- 
efficient of the linear term in B. The divergence and sign change of the magnetic susceptibility 
signal the magnetic instability to the ferromagnetic phase, since its inverse measures the curva- 
ture of the free energy at the origin with respect to the magnetization. Thus quarks near the Fermi 
surface are responsible to the magnetic transition and the spin dependent quark-quark interaction 
and the density of states near the Fermi surface are the key ingredients within FLtQ . 

In this theory, quarks are treated as quasi-particles and thek energy E is regarded as a func- 
tional of the distribution function. FLT is restricted to the low-lying excitations around the Fermi 
surface , where the lifetime of quasi-particles are long enough to be regarded as "particles" with 
properties resembling those of free particles. The single-particle energy is then given by the 
variation of E with respect to the distribution function, incorporating the self-energy; the quasi- 
particle interactions are given by the variation of the single-particle energy with respect to the 
distribution function and they are specified by the set of a few parameters, called the Landau- 
Migdal parameters. It is well-known that there appear infrared (IR) singularities in the Landau- 
Migdal parameters in gauge theories (QED/QCD) since the gauge interaction is infinite range. To 
improve the IR behavior we must take into account the correlation effects for the quasi-particle 
interaction: we must sum up an infinite series of the most divergent diagrams as the self-energy 
term of the gauge field. Thus we take into account the screening effects for the gauge field. 
For example, the Coulomb interaction becomes "finite" range and there is left no singularity in 
the condensed-matter physics The summation of the infinite series or the inclusion of the 
screening effects is also required by the argument of the hard-dense-loop (HDL) resummation 
1 1811 . Consider the particle-hole polarization diagram as a leading-order contribution to the self- 
energy of the gauge field. Since the IR singularities appear for quasi-particles with the collinear 
momenta, the soft gluon should give a dominant contribution. Then the particle-hole polarization 
function n(^) ^ 0{g^jX^) should be the same order of magnitude with the energy-momentum 
squared of gluons ~ {gji)^- Hence we must sum up the infinite series, 1/^^ (l + H/^^ + • •), to 
get the meaningful results. For the longitudinal mode we can see the static screening by the De- 
bye mass and the IR behavior is surely improved. On the other hand, there is no static screening 
for the transverse mode but only the dynamic screening due to the Landau damping lfisll . Thus 
the IR singularities still remain in the Landau-Migdal parameters derived through the exchange 
of the gauge boson of the transverse mode. 

In a recent paper we studied the screening effects for gluons to evaluate the magnetic suscep- 
tibility of quark matter at T = within FLT We have seen that the transverse gluons 
still give logarithmic singularities to the Landau-Migdal parameters, but they cancel each other 
in the magnetic susceptibility to give a finite result. Finally the static screening for the longitu- 
dinal gluons gives the ^'*ln(l /g^) term in the magnetic susceptibility ITs, 21], A similar term is 
also obtained as a correlation effect in electron gas and the ferromagnetic region is diminished. 
However, we also find that this term has an interesting behavior in QCD, depending on the num- 
ber of flavors. Consequently the screening effect does not necessarily work against the magnetic 
instability, which is a different aspect from electron gas. 

In this paper, we extend our previous analysis to the finite-temperature case and figure out the 
non-Fermi-liquid aspect of the magnetic susceptibility. Since thermal gluons should give higher- 



'The high-density effective theory should be also a powerful tool in this context f73 . fT^ . and its relation to FLT has 
been discussed HTH . 
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order terms in T, we only consider the interacting quark system within FLT. We take into account 
the effect of the self-energy of quasi-particles by using the explicit expression given by the per- 
turbation method. We shall see that the longitudinal gluons give little temperature dependence 
due to the static screening, while the transverse gluons exhibit an anomalous temperature depen- 
dence for the magnetic susceptibility through the dynamic screening. Shanker has shown that the 
application of the renormalization group (RG) automatically leads to FLT as a fixed-point, where 
all the quasi-particle interactions are marginal ll22ll . except the attractive superconducting (BCS) 
channel. This argument, however, may be applied to the case of the short-range interaction. Re- 
cent RG arguments have revealed that the quasi-particle interaction through the exchange of the 
transverse gauge bosons is relevant and induces the non-Fermi-liquid effects in gauge theories 
(QED/QCD) 01 H IB [H. Furthermore, it should be important to note that quasi-particle 
interaction is infrared free, so that the effective coupling becomes very weak when the excitation 
energy of the quasi-particles from the Fermi surface approaches to null III6II . The effective cou- 
pling constant is given by the product of the gauge coupling constant and the Fermi velocity 
vp, Qff = g^vp/^lt. Since vp goes to zero on the Fermi surface, the perturbation method should 
be meaningful even if the gauge coupling is not weak. 

The quasi-particle energy et should be given by the root of the transcendental equation, 
£k = £k + ReJ^+(ek,k) with E\i = (|kp + m^)'/^, and the one-loop result of the self-energy 
E+(ek,k) by way of the microscopic many-body technique exhibits an anomalous behavior as 
£k — > M, ReE+(ek) ^ /(Ck — ^)ln(A/|ek — with 7 = g^/9n'^, due to the dynamic screening 
ilA I25I1 . The further analysis of the self-energy by way of the Schwinger-Dyson equation or 
the renormalization group have shown that the anomalous term in the one-loop result is reli- 
able even if 7ln(A/|£k — m|) ^ 1: actually there is no difference between the one-loop result 
and the solution of the Schwinger-Dyson equation 1 16]. Consequently, the renormalization fac- 
tor z+{k) = (1 — dKeL^{(o) / d(o)^^\a,=E^^ vanishes at ^ = kp, which means that the distribution 
function has no discontinuity at e^^^ = /i even at T = 0, or there is no Fermi surface |0, 2^ 271. 



However, the distribution function of the quasi-particles is still the step function at T = 0. Actu- 
ally we have seen that the transverse mode causes no effect for the instability of the Fermi surface 
due to its sharpness 

At r 7^ 0, the Fermi surface is diffused over the width of 0{T) , so that the transverse gluons 
should give rise to an anomalous logarithmic dependence for physical quantities through the 
dynamic screening, as we shall see in the following. For instance, the anomaly in the specific heat 
has been repeatedly discussed in the literatures |24l28, 1, 32, 33]; the transverse gauge 

field gives a T In T term as a leading-order contribution. Analogous non-Ferm i-liq uid behavior 
can be also seen in the context of color superconductivity at high densities 11341 135 , SEt']. We 
shall see that there appears a T^\\\T term in the magnetic susceptibility as another non-Fermi- 
liquid effect i38ll . This temperature dependence is different from that of the specific heat, while 
one may naively expect the same one by considering the density of states on the Fermi surface: 
the leading-order contribution is In T, but it is canceled by another In T term coming from the 
spin-dependent Landau-Migdal parameter to leave the next-to-leading-order term of In T in 
the magnetic susceptibility. 

To extract the temperature dependence correctly we must consider the temperature depen- 
dence of the chemical potential as well. We derive it within FLT and show that the result is 
the same as that given by a microscopic calculation [32]. Properly taking into account the tem- 
perature dependence of the chemical potential, which also includes the In T term, we finally 
present the magnetic susceptibility at finite temperature. 

If the magnetic susceptibility diverges at some density and some temperature, it signals the 

3 



magnetic instability of quark matter to the ferromagnetic phase. Thus we can extract the critical 
densities and temperatures to draw a critical curve on the density-temperature plane. Using our 
results, we shall demonstrate the magnetic phase diagram, which might be useful in considering 
some phenomenological impUcations of magnetism of quark matter. We get several tens MeV 
as the maximum critical (Curie) temperature. This value is similar to that expected during super- 
nova explosion, so that it would be interesting to consider the thermal evolution of magnetars in 
the initial cooling stage. 

In this work we sometimes need the explicit result of a microscopic calculation rather than 
general framework of FLT. In each step of the calculation we try to figure out the difference and 
similarity, compared with the usual discussion of FLT. 

In §2 a framework is presented within Fermi-liquid theory to study the magnetic properties of 
quark matter. In §3 non-Fermi-liquid effects are extracted in the magnetic susceptibility at finite 
temperature. Some specific temperature dependence of chemical potential is also discussed. 
In §4 a magnetic phase diagram is presented on the density-temperature plane, where we can 
also see an importance of the non-Fermi-liquid effect. Some relations between the microscopic 
calculation and the FLT are summarized in Appendix B. 



2. Fermi-liquid theory for quark matter 

Within the Landau Fermi-liquid theory (FLT) we assume a one-to-one correspondence be- 
tween the states of the free Fermi gas and those of the interacting system. Quarks are treated as 
quasi-particles carrying the same quantum numbers of the free particles, and the quasi-particle 
distribution function is simply given by the Fermi-Dirac one, 

«(k,C)= [l+exp(/3(ek,^-M)r' (1) 
with the quasi-particle energy ^ specified by the momentum k and a spin quantum number 

c=±i. 

Since "spin" is not a good quantum number in relativistic theories, we specify the polarization 
of each quark by introducing the spin vector a^, which is a covariant generalization of the spin 
direction ^ (we take it along z— axis) in the non-relativistic case. Actually spin vector is a 
space-like four vector with the constraints, o'^A;^ = and = — 1, and it is reduced to the three 
vector C in the rest frame. It is not uniquely given, but there are many choices for the explicit 
form of spin vector We must choose the most relevant one for the description of ferromagnetic 
phaselQ]. Here we assume the standard form ll39l.l4oll . 

m[ti^ + m) m 

with £i = (k^+m^) 1/2. 

We can define an eigenstate m'?' of the operator, W ■ a, where W'' is the Pauli-Lubanski 
vector, = - \ /^e^vpak" aP" : 

(3) 

m 2 

with ^ — ±L Accordingly we can define polarization density matrix p{k, C,), 

p{k,Q = :^{jt + m)P{a), (4) 
Zm 

by the projection operator, P(a) = l/2-(l + Ys^). 
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2.1. Quasi-particle interaction 

Since the color-singlet magnetization is not directly related to the color degree of freedom, 
we, hereafter, only consider the color-symmetric interaction among quasi-particles that can be 
written as the sum of two parts, the spin symmetric (/^ ^) and anti-symmetric (/^'q) terms; 

/kf ,qC' = /k,q + CC7k,q- (5) 

Since quark matter is color singlet as a whole, the Fock exchange interaction gives a leading con- 
tribution [41]. We, hereafter, take a perturbation technique and consider the one-gluon-exchange 
interaction (OGE) with relatively large coupling constant. One may wonder the perturbation 
method should break in that case. However, the real expansion parameter is not the gauge cou- 
pling constant in the present case; RG argument has shown that the quasi-particle interaction is 
infrared free and the effective coupling becomes very weak for the low-lying excitations around 
the Fermi surface flB']. The effective coupling constant is then given by the product of the gauge 
coupling constant and the Fermi velocity Vf, C^s — g^Vf/4-K, and Vf becomes vanished on 
the Fermi surface. Hence the perturbation method is meaningful even if the gauge coupling is 
not weak. 

For a pair with color index {a,b), the Fock exchange interaction gives a factor {Xa)ah{^a)ba — 
1/2—1/ {2Nc)5ab, which is always positive for any pair Hence the situation is very similar to 
electron gas in QED. Since we are interested in the electromagnetic properties of quark matter, 
only the color symmetric interaction is relevant, which is written as 

I -r^ mm 

/kf ,qC' = 7^ L/kC«,q?'fe = TT TrA^k?,q?' , (6) 

with the invariant matrix element, 

MkCqC' = Vptr(A„/2A„/2)M''^(^,C;^,C'pMv(^-?), (7) 

where M''^ (A;, 1^;^,^') —Xi[y^p{k,C,)'Y^p{q,ll,')]. The matrix elements are easily evaluated with 
the definition of the density matrix p{k^ C,) (Eq. (4)) and summarized in Appendix A. Since the 
OGE interaction is a long-range force and we consider the small energy-momentum transfer 
between quasi-particles, we must treat the gluon propagator by taking into account HDL resum- 
mation [18.1 . Thus we take into account the screening effect, 

D^,{k-q)= P'^.Dr (p) + Pl,D, {p}-^ ^ (8) 

with p = k — q, where [p) ~ [p^ — !!,(;•) )^^, and the last term represents the gauge depen- 
dence with a parameter ^ . P'fj!) is the projection operator onto the transverse (longitudinal) mode, 

P'flV = — gfiV + ^-^-f P^jiv (9) 
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The self-energies for the transverse and longitudinal gluons are given as 



n/(po,p) 
nr(po,p) 



f=u,d.s 



D,f PO 



2uFj IpI 



'I 



(10) 



f=u,d,s 



in the limit po/|p| — *■ 0, with Uf j = kf f/Ef f and the Debye mass for each flavor, mjj ^ = 

g^ji/kfj/ln^ iflill • Thus the longitudinal gluons are statically screened to have the Debye 
mass, while the transverse gluons are dynamically screened by the Landau damping, in the limit 
Po/IpI ^ 0. Accordingly, the screening effect for the transverse gluons is ineffective at T = 0, 
where soft gluons (po/|p| ^ 0) contribute. At finite temperature, gluons with po ^ 0{T) can 
contribute due to the diffuseness of the Fermi surface and the transverse gluons are effectively 
screened. 

2.2. Magnetic susceptibility 

We consider the linear response of the normal(unpolarized) quark matter by applying a small 
magnetic field B. The magnetic susceptibility is defined as 



Xm — ^ Xm ^ H 

f=ii.d.s f=u^d.s 



d(M)f 



dB 



(11) 



B=0 



with the magnetization (M) / for each flavor, where we take B//z. We consider here the spin 
susceptibility, so that the magnetization is given by (M)/ = V^^{J (fixq/L^qf) with /x^ being 
the Dirac magnetic moment. Hereafter, we shall concentrate on one flavor and omit the flavor 
indices because the magnetic susceptibihty is given by the sum of the contribution from each 
flavors. Magnetic susceptibility is then written in terms of the quasi-particle interaction lf5i l20ll . 



Xm 



\ 2 J \+N{T)r 



where go = j\-k\=kF dH^/A-ngoi^) 1^ the gyromagnetic ratio , 

^D(k)C = 2tr[E,p(^,C)]=2(fl, 
and can be explicitly written as 



-flO 



1 



Ek{Ek + m) 



(12) 



(13) 



(14) 



using the spin vector (Eq. (2)). 

A^(r) is the effective density of states on the Fermi surface, and defined by 



N{T) = - 



dn{(0) 
dco^—v{co), 

OCO 



(15) 



^The Debye mass is given as e^p?-Uf / for electron gas in QED. 



where n{co) is the Fermi-Dirac distribution function and v{co) is the density of states of the 
quasi-particles at energy (O, 



v(«) 



2Nc 



I 



Im [^f(«,k)], 



(27r)3 



(16) 



n 



in terms of the retarded Green function of the quasi-particles, J^:f^(a),k) = (o) — Ek — irj) 



Since n{co) d{pL — O)) as T ^ 0, N{T) is reduced to the effective density of states at the 
Fermi surface of the quasi-particles at T = 0, A^(0) = Nck^ / n^Vf with the Fermi velocity, vp — 
{d e]i / d k) k^iij, . Note that —dn{co)/dco is sharply peaked at o = /X for <C 1 which we are 
concerned with, and we can see that only the quasi-particles with the energy Ek — M still gives a 
dominant contribution . 

/" is a Landau-Migdal parameter defined by 



where k^ is defined by e^^ — n and coincides with the usual Fermi momentum kp 3XT — 0. 

3. Non-Fermi-liquid behavior at finite temperature 

We consider the magnetic susceptibility at low temperature. In the usual FLT it should have 
little temperature dependence, while we shall see that the anomalous In T term appears due to 
the dynamic screening effect; the anomalous term is a leading-order contribution at low temper- 
ature. To show such term we carefully analyze the self-energy E+(ek), the derivative of which 
logarithmically diverges as gk — > /i. 

3.1. Average of the density of states near = ^ 

In this subsection, we study the average of the density of state given by Eq. dTSl ). As we 
shall see, we can not use the usual low-temperature expansion because of the singularity of the 
quasi-particle energy at the Fermi surface. Therefore, more careful treatment is needed. 

First, we substitute Eq. ( fTSI l into Eq. ( fTSl ) and change the variable of integration in Eq. (fTsT i 
from ^ to ek(= co); 





(18) 



after angle-averaged over the Fermi surface. 




(19) 




(20) 



'We shall see that the critical (Curie) temperature Tc is actually small compared with mass or Fermi momentum such 
that Tc/moiTc/kf < 1. 
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and 



with Eq = e|k|=o- The quasi-particle energy co satisfies the relation, 

a) = £i,(«)+ReE+(a),^(a))), (21) 

where we discard the imaginary part within the quasi-particle approximation. 

The one-loop self-energy is almost independent of the momentum, and can be written as 

ill El 

ReI+(«,^) ~ ReE+(M,^f ) -^{^{co- fi) In 

\2n^ \(£i — pL\ 

+ ^''\co~^) (22) 

lmL+{(0,k)^^\(0-pL\+0{{(0-pLf) (23) 

around (O ^ fJ. with Cf = {N^ — l)/{2Nc) and Uf = kp /Ej^^. A is a cut-off factor and should 
be an order of the Debye mass. (In this paper, we take A = Mo = \jY.fm^^.) Note that the 
anomalous term in Eq. (|22] | appears from the dynamic screening of the transverse gluons, and 
the contribution by the longitudinal gluons is summarized in the regular function A"^s(a) — /i) 
of 0{g^). Within the approximation given by Eqs. (|22] | and (l23l l. the self-energy is independent 
of spatial momentum k and thus we omit the argument k hereafter The renormalization factor 
z+{k) is then given by the equation, z+{k) = {I — (9ReE+(a))/(9a)|e,=£^)^\ and we have 

z+{k)-' ^ -^f^\n\e^- (24) 

It exhibits a logarithmic divergence as et jJ., which causes non-Fermi-liquid behavior iflill . 
Differentiating Eq. (|2T]) with respect to k, we find 

dk{co) _ <9ReE+(«)^^ £i(<«) ^25) 



d(o \ dm ) k{a)) ' 

Here, we have used (9ReE+ /dk ~ for the self-energy Eq. ( |22] |. 
Eventually, A^(r) is written as. 

As is discussed in Appendix B, we can separate the contribution by the longitudinal gluons Ni{T) 
from N{T). Since the longitudinal gluon exchange is short-ranged by the Debye screening mass, 
it becomes almost temperature independent, 

mT) = '^^fi„ ill) 

with the Landau-Migdal parameter f^. j , 

fh = - ^^Lif [-^4+24] [(1 + ^M^) - 1] , (28) 



where K — Y./ntof/'^kp and 



du 



2 7-1 1 



1 /2 

2 V K- 



(29) 



To evaluate the transverse contribution, Nt{T) ^ N{T) — Ni{T), we only use the transverse 
part in Eq. (l22l i: substituting Eq. (l22l) into Eq. ( |2^ . we obtain the leading-order contribution^ , 



Nr{T) 



(2kj 



CfgW (24-m2) ^2^^/A 
24 A:i "Ir 



^Inf^ 
127r2 V T 



0{g^T^), (30) 



after some manipulation (Appendix C). Nt{T) has a term proportional to InT, which gives a 
singularity at T = 0. This singularity corresponds to the logarithmic divergence of the Landau- 
Migdal parameter /j at T = (Appendix B). Inverting this for T ^ 0, we find 



n-\t) 



{2kj 



fe4 



Cfg^UF (24 -m2). 



72 



4 



-Nl{T) + 0{gV). (31) 



Note that this formula is meaningful only at nonzero temperature; a renormalization group argu- 
ment tells that theory is infrared free in this case and the perturbation analysis of the low energy 
behavior is reliable |16|. We have kept the next-to-leading term {T^lnT term) as well as the 
leading-order term (In T term), since we shall see that the In T term is canceled out by another 
term appearing in the spin-dependent Landau-Migdal parameter /" in the magnetic susceptibil- 
ity. 

3.2. Cancellation of the singular terms at T = 

At finite temperature, the magnetic susceptibiUty is given by 



-1 



(32) 



where ff and ff denote the longitudinal and transverse parts of /" respectively, has the 
logarithmic singularity, which comes from the OGE interaction between quarks with collinear 
momenta. The longitudinal component f" has no IR singularity because of the static screening 
and is almost temperature independent as fj.^ in Eq. ( |27] |. 

The leading-order contribution at finite temperature comes from the transverse component 
f"; it has a logarithmic singularity at T = due to the dynamic screening effect. In this section. 



''We discard here the temperature independent term of O(g^), which cannot be given only by Eq. 1221 . However, 
can recover it by talcing the 7^0 limit later (.see subsection 3.3). 
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we shall see that the logarithmic divergences of ' {T) and //' at T = cancel out each other to 
give a finite contribution to the magnetic susceptibility, f" is given by 



with 



fv.k.ks 



ff = -2NcN-\T) 
' dQ.]i f dD.q 



4n J 4n E.E]^ 



CfN-'g^M''"D,{k-q) 



(33) 



(34) 



where M"" is the spin-dependent component of M" in Eq.(7) (see Appendix A) and we have 
defined Es by Es = E^^^^j^^. It is the dynamic screening part in the propagator D, that gives the 
InT dependence to ff. Therefore, we can put |k| = in the other parts of the integrand in Eq. 
Then, 

/t? r , ,1 

1 + cos + cos 0k COS 0q — COS 9^ — COS 0q . (35) 



M" 



Am? 



The real part of the transverse propagator is 

RtDt{k-q) 



{k~qf 



\l\=ks 



{{k-qfY 
1 



(£k-£.) 

-WW 



H=ks 



- 2^2(1 -cos 0^)3 +c3(£k-£.v)2 

while the imaginary part is given by 



(36) 



(37) 



It is easily shown that the imaginary part gives f only higher-order terms with respect to tem- 
perature and thus we neglect it here. Performing the angular integrals in Eq. ( |34] |. 



48£2 



f" 

Jt;k,ks 



10 + 21n 



cHE^-E, 



c3(£k-£.)2 + 8 

5c(£k-£.v)'/\ ^^-i / 4-c(£k-£.v)'/^ ' 
6 ™ I V3c(£k-£s)2/3 



5;r 
6V3 



c(£'k 



n2/3 



For small \Ek — Es\ ~ it reads 



-1„2 



8£2 



In 



'c{Ek-E,)y^' 



(38) 



(39) 



10 



The integral in Eq. (|33] | can be performed as in the case of N{T), Eq. (20). Changing the variable 
of integration from k to co, we find 



d(ok{(o)E,,(a) 



5ReE+(c(j) 



ft" 



k(a>),ks 



(40) 



Note that dReL+[(o) / d(0 in the integrand does not contribute as long as we calculate to 0{g^) 
as , itself is 0{g^). Integrating Eq. (|40] | as in Appendix C and using Eq. (l3Tl l. we find a 
leading-order contribution at J ^ 0, 



f?-N-\T) 
Cfg^ 



Cfg^ 



1 



7Z^ {2k] 



k^ 



InT 



2t2\ 



0{g'T 



■\nT 



(41) 



Compare Eq. dHJ with Eq. (EB- Since £, = Ep + 0{T^) and k 
shall see, the In T terms cancel each other in the magnetic susceptibility 

n^/{Nck,ii)-N-\T), ff, and N-\T) +/," at fc^ = 1-0 fm"' near T = 



kp + 0{T^) as we 
. In Fig.l, we plot 
As is shown, both 

A^,^ ' (r) and f" are singular at T = 0, but a cancellation occurs so that the sum of the two becomes 
finite even at T = 0. Thus we can take the limit T ^ in the magnetic susceptibility. 

It is also worthwhile to note that Nt{T)+ f" seems to be almost temperature-independent in 
the unit of n^/ (NckslJ.) = 1 in Fig. 1. However, the chemical potential has the implicit tempera- 
ture dependence. We shall discuss this T-dependence in the following. 



Nr\T)+f," 



0.2 



fcF=1 -0 [fm' 



30 40 

r[MeV] 



Figure 1: Cancellation of the logarithmic terms. Tt^ / {Nch^l) ~ N^"^ {T), ff, and N^^(T)+ f" are plotted for A = Mp 
and = 1.0 fm^' in the unit of Tl^ / {N^ks^) = I. f" and N^' (T) are singular at T = 0, but they cancel each other and 
the magnetic susceptibility remains finite even at T = 0. 



3.3. Magnetic susceptibility at T=0 

Before presenting the full expression of the temperature dependence, it would be interesting 
to see how to recover the result at T = |T3 20j|]. We have seen that the logarithmic terms 
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in Nt{T)^^ and ff cancel each other to leave the regular terms in the magnetic susceptibility 
at r = 0. One can easily check that there arises no more singular terms from the longitudinal 
contribution f" atT — 0, since the longitudinal gluons are screened by the Debye mass Mjj- The 
Landau-Migdal parameter f" is given by 



with 



ff = -2NcN-\T) 



{2nr d£k 



fu 



f'l 



I 



CfN;^g^M^^"D,{k~q) 



An J 4n EjEk 



\(l\=l<s 



(42) 



(43) 



in the same way as where m""" is the spin-dependent component of M"" in Eq.(7) (see 
Appendix A). Now the gluon propagator D/ [k — q) has no infrared singularity, so that 



J I —Jl;k,.k, 



d£l^ 



0{g'T 



--H=kF 



(44) 



Consequently /" has little temperature dependence, and can be written as 



ff = - 



2Ejlo{K) + -K{EF-m)\{K) 



K kjloiK) + - ((1 + kMk) - 1) (£f - my \-k'p--{EF- m 



(45) 



with K = L/w^D f/'^kjp at r = l^j- Thus there is no singularity at T = in the magnetic 
susceptibility ( |52] |. 

The magnetic susceptibility at T = is already given in ref. 2oll . 



NckplM) 



= 1 



Un^EUp 



m{2Ef + m) --{Ep+ AEpm - 



2m-^] K-ln- 



(46) 



where Jpauii is the Pauli paramagnetism, ;i;pauii = gf)^^Nckf^/4n , and /io = £'f + 0{g ) is the 
chemical potential at T 0. Using Eqs. ( l3Tl l. (|4TI) . i45[ for the magnetic susceptibility ( |32| | 
and comparing it with Eq. ( |46] | at T = 0, we find that they are different from each other by the 
temperature-independent term Cfg^kf /24-7i^Ef, which is attributed to the contribution by the 
transverse gluons. 

Note that the screening effect for the longitudinal gluons gives rise to the contribution of 
g'*ln(l /g^) ijTs, 21], Obviously this expression ( l46l l is reduced to the simple OGE case without 
screening in the limit K" ^ 0; one can see that the interaction among massless quarks gives a null 
contribution for the magnetic transition. 
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Figure 2: Magnetic susceptibility at T = 0. Tlie solid curve shows the result using simple OGE, while the dashed and 
dash-dotted ones show the screening effects with Nf = l(only s quarks) and with Nf = 3(u, d, and s quarks) respectively. 



In Fig. 2, we plot the magnetic susceptibility at T = O O 19112011 . We take the QCD coupling 
constant as Oo = /An = 2.2 and the strange quark mass rris = 300MeV, inferred from the MIT 
bag modelf42j. We consider here the MIT bag model as an effective model in our energy scale, 
which succeeds in reproducing the low-lying hadron spectra. The coupling constant looks rather 
large, but this value is required for the color magnetic interaction to explain the mass splitting of 
hadrons with different spins; e.g. for nucleon and A isobar. We think this feature is relevant in our 
study, because the coupling constant is closely related to the strength of the spin-spin interaction 
between quarks in this model. Moreover, the quark density in the MIT bag model is moderate, 
^ 0.25fm^^, which is the similar one we are here interested in. Note that the perturbation method 
should be still meaningful even for this rather large coupling, since the renormalization-group 
analysis has shown that the relevant expansion parameter is not the gauge coupling constant 
but the product of with the Fermi velocity vp, which always goes to zero as one approaches 
to the Fermi surface iflill . 

One can see that the magnetic susceptibiUty for the simple OGE without screening diverges 
around kp = 1.3 fm This is consistent with the previous result for the energy calculation. 
|fl9 , 20]. One may expect that the screening effect weakens the Fock exchange interaction so 
that the critical density get lower once we take into account the screening effect. However, this 
is not necessarily the case in QCD. The screening effect behaves in different ways depending on 
the number of flavors. Compare the results for the Nf = 3 with the one for Nf = 1. In the case 
of Nf = 1 , K" < 2 the screening effect works against the magnetic phase transition as in QED. 
However, for Nf = 3, fc > 2 so that the critical density is increased. Consequently the screening 
effect does not necessarily work against the magnetic instability, which is a different aspect from 
electron gas 1I20I1 . 

3.4. Temperature dependence of the chemical potential 

The chemical potential jX in Eq. OTT l implicitly includes the temperature dependence. To 
extract the proper temperature dependence in Xm we must carefully take into account the tem- 
perature dependence of /i ; here we use the thermodynamic relation = —{dF/dn)\T with the 
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free energy F = E — Ts. 

In the following we only consider the temperature variation 5T, so that the contribution of 
the longitudinal gluons can be well discarded due to the Debye screening. Accordingly it is 
sufficient to take into account the contribution from the transverse gluons in Eq. ( |22] |. We shall 
see that includes In T term due to the dynamic screening effect for the transverse gluons, 
besides the usual term. 

Using the fact that the temperature variation of the free energy equals —s8T, we calculate 
the quasi-particle entropy per unit volume s, 

^ = -^'E /^Kk,C)ln«(k,C) + (l-«(k,C))ln(l-«(k,C))] (47) 

for each flavor, and then integrate it with respect to temperature to get F. The temperature 
variation of s is given as 

Ss = -^N,.'£ J ^(e,^-^)5„(k,C) (48) 



through the variation of the distribution function. 



Mk,C)^^ 



(49) 



Since the term involving Se^^ — gives at least InT in the entropy at low temperature (It], 
we shall see that the leading-order contribution is given by the term with the explicit temperature 
variation dT, 

5s=-Ncy / 7 — £kC-M — • (50) 



Using Eq. ( |25] |. 5s is recast into the following form. 



Changing the variable o by y s.t. (O — jX =yT and using Eq. ( 1221 ). we have 
as in Appendix C. Thus the entropy is 

Note that this result is the same as the one given in ref. ll32i] . and the non-Fermi-liquid behavior 
of the specific heat cy = T{ds/dT)v = s + 0{g^T) is also reproduced. The free energy F is then 
given by the integral of s with respect to temperature. 



F^Eo + ^-^T' + ^^^^T'\nl^-]+0{g'T'), (54) 
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where Eq is the ground state energy at T = 0. Through the thermodynamic relation /X = — {dF /dn) 
we finally obtain 



klEp 



In 



2t2\ 



0{g'T 



(55) 



Note that this relation is also given by considering the total number density at finite temperature 



3.5. Temperature dependence of the magnetic susceptibility 

One can easily find that temperature-dependent terms of is the same as those of /i(r) 
within the order we are interested in, noticing that E^ = fi — ReE+ (jj.) and Rei;+ (jJ.) — ReT.+ (/io) + 
0{g^T^). Using ks = \/E} — m^, we find the temperature dependence of ks. 



ks — kf 



2kl - 



k^ 



Now we can rewrite Eqs. dSTT i and ( |4T1 ) in terms of kf and Ef; 



0{g^T^). 



(56) 



n-\t) = 



Nckpfio I 6kj. 
Cfg^Uf 2kp — m^ 



2Ei-m^ + -^ T 



£■2 



12 



kp 



r^in 



Cfg^UF 

^2^'\t 



-NiiT)+0{g^T^) 



(57) 



and 



f" = 



Cfg' 



UNcEffio 



F F 



(58) 



Taking into account this temperature dependence of the chemical potential and using Eqs. (|57] | 
and ( fSSl ), we find the temperature dependence of SXm^ , 



^Pauli 



64 



2E^p~m'- + —y I 



12k'pEl 



{2kt 



-klm^ + m'^)T^\n{ - 



0{g'T') 



(59) 



Finally we obtain the magnetic susceptibility at finite temperature by adding 5;)^^' to Eq. (|46]) . 



iXM/XPauuY 



64 



Cfg^ 



Un^Efkf 

^4 



1 2 1 

m(2EF + m) - -{Ep + AEpm - 2m^)Kln - 



2Ep-m^ + — ]T 
Ep 



Cfg^UF {2kj + klm^ + m^] 



72 



-0{g^T^) 



(60) 



There appears In T dependence in the susceptibility at finite temperature besides the usual 
dependence ll20ll . This corresponds to the TlnT term in the specific heat 1 24 , 28 , 2^ 3^ 31 , 
32|,|33|] and is a novel non-Fermi-liquid effect in the magnetic susceptibility. At low temperature, 

15 



ln(A/r) > so that the T^lnT term gives positive contribution to Xm^ ■ Therefore, both T- 
dependent terms in Eq. (l60t work against the magnetic instability. 

In Fig. 3, we show the temperature dependence of the magnetic susceptibility at kf — 1 .0 fm^ ^ 
using Eqs.(|59b and ( l60l ). The thin solid, dash-dotted, heavy solid, and dashed curves indicate 
SXMiT)/XPauii, Xm{T = 0)IXv^u\u ZM(r)/ZPauii and {kFlla)/{ksPi) for A = Mb respectively. 
The temperature at which the heavy solid curve cross the zero corresponds to the critical(Curie) 
temperature. Temperature dependence of the magnetic susceptibility, i.e. 5xm originates from 
the transverse gluons. {kfpLQ)/{kspL) is rewritten as 

and corresponds to the first term in Eq.dsTli. Once we take into account the temperature de- 
pendence of the chemical potential, this term depends on temperature and adds temperature 
dependence other than the explicit temperature dependence to the magnetic susceptibility. We 
also plot the magnetic susceptibility with the temperature dependence of the chemical potential 
disregarded. If we ignore the temperature dependence of the chemical potential, the magnetic 
susceptibility hardly depend on temperature. The temperature dependence of the magnetic sus- 
ceptibility almost comes from that of the chemical potential. 
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Figure 3: Temperature dependence of the magnetic susceptibility at i:;r = 1.0 fm ' . SXnf\T) /Xp.^^\i, Xm^C^ ~ //Cpauli ' 
Xm' (^)//tpaIiii> {^FPn)/{ksli) are plotted for A = Mo- The dotted curve is the magnetic susceptibility in which the 
temperature dependence of the chemical potential is neglected. The parameters are the same as those used in Fig. 1. 
(T) changes its sign at T = Tc, which gives a critical (Curie) temperature. 



4. Magnetic phase diagram of QCD on the density-temperature plane 

In this section, we show some numerical results for the magnetic susceptibility at finite tem- 
perature. We consider a symmetric quark matter with three flavors in equal populations, massless 
u^d quarks and massive s quarks. We use the same values for m and as in Fig. 2. 



16 




kF [1/fm] 
(b) 

Figure 4: Magnetic susceptibility at finite temperature. The upper panel shows the result for the full expression. Eg. 1601 . 
while the lower one shows the result without the T~lnT term in Eq. )60t . The dotted, dashed, dash-dotted, and solid 
curves show the results at r=0, 30, 50, and 60 MeV respectively. The diamonds indicate the critical density, atT = Q 
and the open (filled) triangles, squares, and circles indicate kpi{k'p2) for T =30, 50, and 60MeV respectively. 
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In Fig.4(a), we plot the magnetic susceptibility given by Eq. i&jj . At T=0, the magnetic 
susceptibility is positive at higher densities and the quark matter is in the paramagnetic phase 
there. At the critical density where the magnetic susceptibility diverges (kp ^ 1.6fm^'), there 
occurs a magnetic phase transition from the paramagnetic phase to the ferromagnetic phase and 
the quark matter remains ferromagnetic below k^p. 

At T =30 MeV, there appear two critical densities at which the magnetic susceptibility di- 
verges. We denote these densities kp^ and kp2 (k^pi < kp^. In this case, kpy ~ 0.4fm^^ and 
k'^P2 ~ 1.5 fm^^ At densities below fc^j and above k'p2, the magnetic susceptibility is positive, 
which corresponds to the paramagnetic phase, on the other hand, at densities between two critical 
densities, it becomes negative corresponding to the ferromagnetic phase. 

At T =50 MeV, there ai-e still two critical densities ~ 0.7fm"' and kp2 ~ 1.3fm"'), but 
the range between these two densities becomes narrower than at T =30 MeV. 

At T =60 MeV, there is no longer divergence in the magnetic susceptibility and the quark 
matter is paramagnetic at any density. 

To figure out the non-Fermi-liquid effect, in Fig. 4(b), we depict the magnetic susceptibility 
without the T'^ In T term in Eq.(l60ll for comparison: only T^- dependence on temperature is taken 
into account. At T = 0, there is no difference from the curve shown in Fig. 4(a) and there is little 
difference between Fig.4(a) and (b) even at T =30 and 50 MeV. However, at T =60MeV, there 
is a region where the spontaneous magnetization is still realized in Fig.4(b), whereas the quark 
matter is paramagnetic at any density in Fig. 4(a). From these considerations, the non-Fermi- 
liquid effect lowers the Curie temperature to some extent and affects the magnetic property of 
quark matter in contrast with that of the anomalous specific heat {oc T In T) in a non-relativistic 
electron gas S El S EH . 

We show a magnetic phase diagram of QCD on the density-temperature plane in Fig. 5. The 
four curves corresponds to the critical curves given by Eq.(l60]l under four different assumptions: 
below the curves the quark matter is in the ferromagnetic phase, while it is in the paramagnetic 
phase above the critical curves. The magnetic transition occurs on the critical curves. 

For the solid curve, we have used the full expression Eq.(l60l). on the other hand, for the 
dashed, dash-dotted, and dotted curves, we have ignored the dynamic screening(/.e. the T^\nT 
term), static screening(/.e. the fclnfc term), and both of the two screenings in Eq. ( l60b respec- 
tively. 

Compare the result with the full expression (|60] | with the one without the non-Fermi-liquid 
effect i.e. T^lnT dependence. In the case without the T^lnT term, the ferromagnetic phase 
can be sustained till over T = 60 MeV, while it can be at most T = 60MeV including In T 
dependence. It turns out that the dynamic screening works against the magnetic instability and 
can reduce the ferromagnetic region in the phase diagram up to a point, but this effect is not so 
large. 

The dash-dotted curve is the result without the static screening or fcln K term in Eq.(l60t. The 
static screening effect works in favor of the magnetic instability to enlarge the ferromagnetic 
region. As discussed in i20tl . it depends on the number of flavors whether the static screening 
works for the ferromagnetism or not, which is peculiar to QCD. 

The maximum Curie temperature T^^^ is around 60MeV, which is achieved at kf ~ I . Ifm^ ' . 
Note that this is still low temperature, since T^^^/kf ^ 1. Thus our low-temperature expansion 
is legitimate over all points on the critical curve. One of interesting phenomenological implica- 
tions may be related to thermal evolution of magnetars; during the supernova expansions tem- 
perature rises up to several tens MeV, which is so that ferromagnetic phase transition may occur 
in the initial cooling stage to produce huge magnetic field. 
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Figure 5 : Magnetic phase diagram in the density-temperature plane. The soHd, dashed, dash-dotted, dotted curves show 
the results for the full expression Eq. )60t , the one without the T^lnT term, without the KlnK term, and without the 
T- In T and KlnK terms in Eq. 160) . The open (filled) circle indicates the Curie temperature 3X kf = 1. 6) fm^ ' while 
the squares show those when we disregard the T-lnT dependence. 



In Fig|6l we plot the contributions to the magnetic susceptibility from the temperature inde- 
pendent, , and T^lnT terms for kf = 1.1 fm^^ and kf = 1.6 fm^'. We represent the and 
T^lnT term in Eq. ( l60b by 5Xj2 ^"d SXj2i„f respectively. Compare the contributions from 
the two temperature-dependent terms. In either of the cases, the term is dominated over the 
r^lnT term except at extremely low temperature. At sufficiently low temperature, the T^lnT 
term should exceed the term, although we cannot see such a situation in Fig. 6. Now, we esti- 
mate the reversal temperature, Trev, at which dXj2 starts to exceed dXj^^j- This temperature is 
easily obtained from Eq. 



: Aexp 



■ Tye\\c exp 



12n^ 
Cfg^uf 

llTt^ 



2kl-m^ + 



24- 



- k^pirP- - 



Cfg^Uf 



(62) 



where Trev.c is the reversal temperature of the specific heat, at which the T term and the TlnT 
term in the specific heat or the entropy, Eq.(l53Tl are equal in magnitude. Trev is quite low for the 
densities we are interested in. (For example, Trev = (9(0.01 )MeV for kf = 1.1 fm^' and Trev — 
0(1) MeV for = 1.6 fm^'.) As is shown in Fig. 6, the T^lnT dependence is too important to 
neglect, although the reversal temperature of the magnetic susceptibility is quite low. 
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Figure 6: Contributions to tlie magnetic susceptibility. The left panel shows the result for if = 1 . 1 fm and the right one 
is for = 1.6 fm"' in the unit of XPmYi = 1- The solid, dashed, dash-dotted, and dotted curves show SXj^^ {T), SXj2 ' 
SXjifj^Y' ^'^^ ^'t^jw'C) respectively. The circles(squares) indicate the critical(Curie) temperatures including(ignoring) 
the In T dependence. SXj2fji rneans the (T~ In T) term in the magnetic susceptibility, Eq. 160) . 



5. Summary and concluding remarks 

We have discussed some magnetic aspects of quark matter by evaluating the magnetic sus- 
ceptibihty at finite temperature based on QCD, where the screening effects for gluons are taken 
into account. Different from the usual FLT, we cannot use the low temperature expansion be- 
cause the quasi-particle energy have an anomalous term on the Fermi surface. Carefully studying 
the temperature effects, we have found that the magnetic properties of QCD exhibit interesting 
features reflecting the gauge interaction: the dynamic screening of the transverse gluon field 
gives an anomalous T^lnT contribution to the magnetic susceptibility [38], while the contribu- 
tion by the longitudinal gluons is almost temperature independent due to the Debye mass. It may 
be interesting to recall that the dynamic screening has no effect but the static screening gives 
the term proportional to M^lnM^' &\.T — 0[20] with the Debye mass Mo\ Mo works as an in- 
frared (IR) cutoff to remove the infrared singularity in the quasi-particle interaction exchanging 
the longitudinal gluons, while the IR singularity still remains in the quasi-particle interaction 
due to the transverse gluons. At finite temperature, the energy transfer of order T is allowed 
among the quasi-particles in the vicinity of the Fermi surface, so that temperature itself plays a 
role of the IR cutoff through the dynamic screening for the quasi-particle interaction. The log- 
arithmic temperature dependence appears in the magnetic susceptibility as a novel non-Fermi- 
liquid effect, and its origin is the same as in the well-known TlnT dependence of the specific 
heat fM. 28, 2^ 3^ 31, 32, 33]. The dynamic screening also gives InT terms for the spin- 



asymmetric Landau-Migdal parameter as well as the density of states on the Fermi surface as the 
leading-order contribution, but they exactly cancel each other in the formula of the magnetic sus- 
ceptibility. Hence In T term becomes the leading-order contribution. The anomalous In T 
term works against the magnetic instability, as well as the usual T~ term. 

We demonstrated how the temperature effects work by drawing the magnetic phase diagram 
on the density-temperature plane. We also figured out the role of the non-Fermi-liquid effect 
in the phase diagram. We have seen that the In T term should have a sizable effect on the 
magnetic phase diagram in the temperature-density plane . 

We have seen some relations between the microscopic many-body calculation and the FLT; 
the Landau-Migdal parameter /j or for the spin symmetric interaction is closely connected 



^It should be worth noting that the anomalous T^lnT term in the magnetic susceptibility could be also observed in 
electron gas. 
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with the quark self-energy and thereby the renormahzation factor z+ikp). It diverges at T = 
due to the transverse mode, but another divergence coming from the spin dependent interac- 
tion cancel it to give a finite result for the magnetic susceptibility. At finite temperature all the 
Landau-Migdal parameters become finite. We have also derived the temperature dependence 
of the chemical potential within FLT to get the same result as that given by the microscopic 
calculation. 

The magnetic susceptibility is a powerful tool to study not only the magnetic properties of 
QCD but also the magnetic transition in quark matter However, it only gives the information 
about the curvature at the origin of the effective potential (free energy) with respect to magne- 
tization. We have discussed the magnetic phase diagram in QCD, assuming that the magnetic 
transition is of the second order or the weakly first order. If the phase transition is of the first- 
order, we must explore the global behavior of the effective potential. The critical density or the 
critical (Curie) temperature then would be larger than that in this paper. It would be an interesting 
possibility to be further investigated l43[ 14411 . 

Here we have only considered the Fock exchange interaction in a perturbation method to 
reveal the interesting features of the magnetic properties within gauge theories. We have used 
a rather large coupling constant to demonstrate the specific features of the magnetic phase tran- 
sition and the magnetic phase diagram, referring to the bag model parameter, where the color 
magnetic interaction gives a spin-spin interaction to split the mass degeneracy, e.g. for nucleon 
and A isobar. To be more realistic, we must take into account the non-perturbative effects, since 
some non-perturbative effects, like instantons |43, 45], are still relevant at low and intermediate 
densities. Actually some part of the mass splitting should be attributed to such non-perturbative 
effects, which in turn affect the magnetic properties of quark matter. Although the coupling con- 
stant should be small in this case, the non-Fermi- liquid behavior still survives as a qualitative 
effect. 

In high-density QCD color superconductivity (CSC) has been extensively studied by many 
authors [46]. We have completely discarded its possibility in this paper to concentrate in the 
magnetic properties. It would be interesting to study the interplay of superconductivity and 
magnetism in quark matter A first attempt in this direction has been already done j^, but more 
studies are needed to this end. 

Finally, it should be worth noting that magnetic properties of QCD or its magnetic instability 
may be related to physics of compact stars, especially magnetars or primordial magnetic 
field in early universe, where one may expect the QCD phase transition [47]. Our phase diagram 
may have a direct relevance to the thermal evolution of magnetars; the evolution path can be 
plotted on the temperature-density plane, and it may traverse the critical line in the initial cooling 
stage just after the supernova explosion. 
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Appendix A 

The matrix elements for two quarks with ^ and ^' 

M«0 = tx[fp{k,Qfp{q,C% 

M'J = tr[Yp{k,QY'piq,C)], (63) 
can be easily evaluated. A straightforward manipulation, using Eq. (3), gives 

= -^{2koqo-k-q + m^+{np--q-k){2aobo-a-b) 

— [—laoqok ■ b + Ikoq^a ■ b — 2koboa -q + k-ba-q]} , (64) 

and 

M'-i = --^{{\~a-b)k^i + {m^ -k-q)^i + a-qk^j +b-k^i 

{{n?--k-q)(\-a-b)-k-bq-a\}, (65) 

with a symbol, k^qi = k'q^ + k'q', where the spin vector is given by Eq. (2) and similarly bfi 
Appendix B 

We present here some relations between the self-energy Z+ and the Landau parameter. The 
Fermi-hquid interaction of OGE is calculated at T = as 

/kf.qC l|k|=|qh^^ = [-M^'^Diico - M,k - q) +M"A(co - M,k - q)] , (66) 

where the matrix elements Af^^ are given in Appendix A, and the propagator renders 

1 



A(a)-M,k-q)^-^^^^^3^^^^, (67) 



for the longitudinal gluons with K" being the parameter in terms of the Debye mass, K = / Ik^ = 
Y,finjjj-/2kp. The propagator for the transverse gluons renders 

A(«-Al,k-q)||k|=|q|=^^ r T, (68) 



2kl 



l-COS0jp + 'c3/^7^-^^^^T72 



with c'^ = -^(T,f '""°'f"'''-f ^ in the soft-gluon Umit, \co-n\/n<^l. 



It is rather easy to see the contributions to the Landau-Migdal parameters from the longitudi- 
nal gluons, since they have no IR singularity. For the spin-averaged interaction, we have 



2 

) m 



fl = -iV-iC/g2 Ij^j^ M«o^A(ek-M,k-q)||k|=|q|=*^> (69) 
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with 



^00. ^ 1 ,,00 = -1, (4 + 4 COS + m-) . (70) 



Accordingly the Landau-Migdal parameter, f^^ reads 

= - ^^if K + 24] [{1 + kMk)-\]. (71) 

Using the relation of the density of state near the Fermi surface and the spin-averaged Landau- 
Migdal parameter ff'l, we find 1 contributes to A^(0) at T = 0. Since there is little tempera- 
ture dependence for the short-range interaction, we have Ni{T) ~ NckfEf /{37:^)fi.^. Recalling 
Eq. ( l22b . we can make sure that the longitudinal contribution is implicitly included in E+ through 

For the transverse gluons, we must carefully calculate the Landau-Migdal parameters, since 
they never receive the static screening but the dynamic one due to the Landau damping. Generally 
Df includes the imaginary part proportional to the squared Debye mass, but it only gives a higher- 
order contribution. The real part of the Fermi-liquid interaction now reads 

/,k?,q?'||k|Hqh., - - ' /z iigl [^"ReA(a)-Ai,k-q)]|,|^|^j^,^ . (72) 

F 

For the spin-dependent Landau parameter we have ah*eady seen in §3 that 
as a leading-order contribution. For the spin averaged interaction, 

2 

// = -A^-iC/g2|_ jij^ M'"ReA(ek-M,k-q)||k|=|q|=^^, (74) 

with 

M'" = -y M" = 

4 2m2 

One can see that the relevant Landau parameter f^.y also contains InjCk — /x| due to the angular 
integral, f^.^ is defined by 

/;.i = ^|rf(cos0g^)cos0g^/;^ (76) 



//' ~ ' lim Inlfik - M I (73) 



l+-(l-COS0j;,) 



(75) 
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Consider the integral, 



d{cos cos e^ReDfM'" 
1 /■2 

- / (l-x)(2+x)3— , 

OT^ Jo + C-^ 



1 

8^ 



4 + |ln 
o 



(2 + c 



^ tan 



4-c ;r 



V3 V V3c 6 



— In 
6 







(2 + c)3 





tan 



-1 



4-c TT 



V3 V \/3c 6 



(77) 



with c = c(ek — /i)^/^. For a small — it reads 



1 

8^ 



-21n 



Thus the leading-order contribution to f^.^ is eventually given by 



ft- 1 = lim /r;ik, 



(78) 



(79) 



4El 



ln|ek-;u| or 



N{o)e 



Cfg^kp 



lim ln|£k — ^1 



up to 0{g^). Comparing Eq. dSOl l with Eq. (|22] |. we can see 



FA = Hm '-^^ 



(80) 



(81) 



in the leading order. 

Consider the renormalization factor on the Fermi surface, which measures a discontinuity of 
the momentum distribution of the bare particles |0, 2^ 27]. We find 



-1 



0. 



The Fermi velocity is given by 



Moreover, one can see that 



-1 _ ^ 



VF = —z+ikp). 



(82) 



(83) 



(84) 



In the non-relativistic version, this relation is nothing but the relation between the effective mass 
m* and the renormalization factor z+ {kp), m* — z^^{kf). 
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At r 7^ 0, the contribution of /,![ to N{T) can be estimated by evaluating 

d^k dn{ek) .f,\ik _ Cfg^UF 



■2Nc 



(2;r)3 dsk 3 



12;r2 



■ inr, 



(85) 



which is reduced to 1 /3F^''.^ at T = 0. Thus we can see that the In T term in Eq. dSOb is just given 
by the Landau-Migdal parameter f^.^ . 

Appendix C 

Here we derive Eq. ( l30l ). Substituting Eq. (l22l l into Eq. ( |26] | and changing the variable CO by 
the dimensionless one y through co — jJ. — yT, we have 



N{T) ^ § dyk{^ +yT)E,^^+..r) ( 1 - '^'^^^"'^ 

7l J —oo \ 



dco 



m=n+yT , 



(e.v + l)2' 



(86) 



where we put (Eq — /x)/r ^ — oo in the lower limit of the integral. Now we expand the each term 
in the integrand to find the leading-order contribution at low temperature. We thus have 



1 



<9ReE+(a)) 



d(i) 



a=n+yT 



1-^infz: 



and 



(87) 



(88) 



in the integrand, where the odd-power terms of y is discarded since they never contribute to the 
integral with respect to y. Finally we have 



+■ 



2ki 



(e>' + l)2- 



(89) 



The integral in Eq. ( |89] l is easy to give Eq. (l30b . 



References 



[1] O. Philipsen, arXiv:0710.1217 

[2] P. Braun-Munzinger and J. Stachel, Nature 448 (2007) 302. 
[3] For recent reviews, 

F. Weber, Prog, in Part, and Nucl. Phys. 54 (2005) 193. 

D. Page and S. Reddy, Ann. Rev. of Nucl. and Part. Sci. 56 (2006) 327. 

J.M. Lattimer and M. Prakash, Phys. Rep. 442 (2007) 109. 
[4] T. Schaefer, arXiv:hep-ph/0509068 

P. Braun-Munzinger and J. Wambach, arXiv:0801.4256' 
[5] G. Baym and C.J. Pethick, Landau Fermi-Liquid Theory (WILEY- VCH, 2004). 

A. A. Abrikosov, L.P. Gorkov and I.E. Dzaloshinski, Methods of Quantum Field Theory in Statistical Physics 

(Prentice-Hall. Inc., 1963). 

A.B. Migdal, Theory of finite Fermi systems (Intersci. Pub., 1967). 

25 



[6] P. Nozieres, Theory of Interacting Fermi Systems (Westview Press, 1997). 

D. Pines and P. Nozieres, The Theory of Quantum Liquids (Perseus books Pub., 1999). 
[7] G. Baym and S.A. Chin, Nucl. Phys. A262 (1976) 527. 
[8] T. Tatsumi, Phys. lett. B489 (2000) 280. 

T. Tatsumi, E. Nakano and K. Nawa, Dark Matter, p. 39 (Nova Science Pub., New York, 2006). 
[9] E. Nakano, T. Marayama and T. Tatsumi, Phys. Rev. D68 (2003) 105001. 

T. Tatsumi, E. Nakano and T. Maruyama, Prog. Theor. Phys. Suppl. 153 (2004) 190. 

T. Tatsumi, T. Maruyama and E. Nakano, Superdense QCD Matter and Compact Stars, p. 241 (Springer, 2006). 
[10] A. Niegawa, Prog. Theor. Phys. 113 (2005) 581. 

[11] P.M. Woods and C. Thompson, Soft gamma ray repeaters and anomalous X-ray pulsars:magneter candidates. 
Compact stellar X-ray sources, 2006, 547. 

A. K. Harding and D. Lai, Rep. Prog. Phys. 69 (2006) 2631. 
[12] F.Bloch,Z. Phys. 57(1929) 545; 

[13] C. Herring, Exchange Interactions among Itinerant Electrons: Magnetism A' (Academic press. New York, 1966) 

K. Yoshida, Theory of magnetism (Springer, Berlin, 1998). 
[14] T. Moriya, Spin Fluctuations in Itinerant Electron Magnetism (Springer, 1985) 
[15] D.K. Hong, Phys. Lett. B473 (2000) 118; Nucl. Phys. B582 (2000) 451. 
[16] T. Schafer and K. Schwenzer, Phys. Rev. D70 (2004) 054007; 1 14037. 
[17] S. Hands, Phys. Rev. D69 (2004) 014020. 

[18] J.L Kapusta, Finite-temperature field theory (Cambridge U. Press, 1993). 

M. Le Bellac, Thermal Field Theory (Cambridge U. Press, 1996). 
[19] T. Tatsumi, Exotic States of Nuclear Matter (World Sci., 2008) 272. 
[20] T. Tatsumi and K. Sato, Phys. Lett. B663 (2008) 322. 
[21] K.A. Brueckner and K. Sawada, Phys. Rev. 112 (1957) 328. 

B. S. Shastray, Phys. Rev. Lett. 38 (1977) 449. 
[22] R. Shanker, Rev. Mod. Phys. 66 (1994) 129. 

[23] C. Nayak and F. Wilczek, Nucl. Phys. B430 (1994) 534. 

J. Polchinski, Nucl. Phys. B422 (1994) 617. 
[24] D. Boyanovsky and H.J. de Vega, Phys. Rev. D63 (2001) 034016; 1 14028. 
[25] C. Manuel and Le Bellac, Phys. Rev. D55 (1997) 3215. 

C. Manuel, Phys. Rev. D62 (2000) 076009. 

[26] J.M. Luttinger and J.C. Ward, Phys. Rev. 118 (1960) 1417. 

J.M. Luttinger, Phys. Rev. 119 (1960) 1153. 
[27] J.M. Luttinger, Phys. Rev. 121 (1961) 942. 

[28] T. Holstein, R.E. Norton and P Pincus, Phys. Rev. B8 (1973) 2649. 
[29] M.Yu. Reizer, Phys. Rev. B40 (1989) 11572; B44 (1991) 5476. 
[30] J. Gan and E. Wong, Phys. Rev. Lett. 71 (1993) 4226. 

[31] S. Chakravarty, R.E. Norton and O.R Syljuasen, Phys. Rev. Lett. 74 (1995) 1423. 
[32] A. Ipp, A. Gerhold and A. Rebhan, Phys. Rev. D69 (2004) 01 1901. 

A. Gerhold, A. Ipp and A. Rebhan, Phys. Rev, D70 (2004) 105015; PoS (JHW2005) 013. 
[33] A. Gerhold and A. Rebhan, Phys. Rev. D71 (2005) 085010. 
[34] D.T Son, Phys. Rev. D59 (1999) 094019. 
[35] T. Schafer and E Wilczek, Phys. Rev. D60 (1999) 1 14033. 

[36] W.E. Brown, J.T Liu and H. Ren, Phys. Rev. D61 (2000) 1 14012; D62(2000) 054013. 
[37] Q. Wang and D.H. Rischke, Phys. Rev. D65 (2002) 054005. 
[38] T. Tatsumi and K. Sato, Phys. Lett. B672 (2009) 132. 

[39] C. Itzykson and J.-B. Zuber, Quantum Field Theory (McGraw-Hill Inc., 1980). 
[40] V.B. Berestetsii, E.M. Lifshitz and L.P. Pitaevsii, 

Relativistic Quantum Theory (Pergamon Press, 1971). 
[41] K. Ohnishi, M. Oka and S. Yasui, Phys. Rev. D76 (2007) 097501. 
[42] T.A. DeGrand, R.L. Jaile, K. Johnson and J. Kiskis, Phys. Rev. D12 (1975) 2060. 
[43] M. Inui, H. Kohyama and A. Niegawa, arXiv:0709.2204 
[44] K. Pal, S. Biswas and A.K. Dutt-Mazumder, arXiv:0809.0404l 

[45] S.-il Nam, H.-Y. Ryu, M.M. Musakhanov and H.-C. Kim, arXiv:0804.0056' 

[46] For a recent review, M.G. Alford, K. Rajagopal, T. Schaefer, A. Schmit, a rXiv:0709.4635l 
[47] D. Boyanovsky, H.J. de Vega, D.J. Schwarz, Ann. Rev. Nucl. Part. Sci. (2006) 441. 



26 



